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ABSTRACT 

We present a survey of the photometric and kinematic properties of 39 nearby, nearly face-on 
disk galaxies. Our approach exploits echelle-resolution integral-field spectroscopy of the Ha regions, 
obtained with DensePak on the WIYN 3.5m telescope Bench Spectrograph. This data is complemented 
by H I line-profiles observed with the Nangay radio telescope for 25 of these sample galaxies. Twelve 
additional line-widths are available for sample galaxies from the literature. In this paper, we introduce 
the goals of this survey, define the sample selection algorithm, and amass the integral field spectroscopic 
data and Hi line-widths. We establish spatially-integrated Ha line-widths for the sample. We test 
the veracity of these spatially-integrated line profiles by convolving narrow-band imaging data with 
velocity field information for one of the sample galaxies, PGC 38268, and also by comparing to Hi line 
profiles. We find H I and Ha line profiles to be similar in width but different in shape, indicating we 
are observing different spatial distributions of ionized and neutral gas in largely axisymmetric systems 
with flat outer rotation-curves. We also find vertical velocity dispersions of the ionized disk gas within 
several disk scale-lengths have a median value of 18 km s _1 and an 80% range of 12-26 km s _1 . This is 
only a factor of ~2 larger than what is observed for neutral atomic and molecular gas. With standard 
assumptions for intrinsic and thermal broadening for Ha, this translates into a factor of three range 
in turbulent velocities, between 8 and 25 km s _1 . 

Subject headings: galaxies: structure; galaxies: kinematics; galaxies: spiral 



1. INTRODUCTION 

Two fundamental and intertwined cosmological ques- 
tions are how galaxies form and come to have their 
present structure. The current approach to answering 
these questions involves comparing numerical and semi- 
analytical structure-formation models to a suite of crit- 
ical observations of low and high redshift galaxies (e.g., 
Navarro, Frenk, & White 1996, 1997; van den Bosch 
2000). For nearby disk galaxies one of the salient obser- 
vations which tests and constrains structure-formation 
models is a measurement of halo shape. 

An increasing number of cosmological simulations show 
that dark matter halos are naturally triaxial (Davis et al. 
1985; Frenk et al. 1988; Dubinski & Carlberg 1991; 
Jing et al. 1995; Jing & Suto 2002), but other physi- 
cal processes, such as the dissipative infall of gas dur- 
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ing formation, may transform the halo into an oblate 
spheroid (Dubinski 1994). If dark matter halos are non- 
spherical, then this should have an observational impact 
on the luminous matter in disks. Theory predicts a num- 
ber of observational consequences including effects on (i) 
the orbits in lopsided galaxies (Levine & Sparke 1998; 
Jog 1999; Noordermeer, Sparke, Levine 2001; Jog 2002); 
(ii) the formation and structure of bars and oval distor- 
tions (Athanassoula 2003); and (iii) non-circular motions 
in dark matter dominated low surface brightness (LSB) 
galaxies (Hayashi et al. 2004) . Triaxial halos also should 
beget elliptical disks (Franx & de Zeeuw 1992, Jog 2000). 
Because disk asymmetries may be linked to the shape 
of the halo, it is interesting to measure deviations from 
axisymmetry. Indeed, many galaxies appear to have lop- 
sided light distributions (e.g., Zaritsky & Rix 1997; Korn- 
reich, Haynes, Lovelace 1998; Conselice, Bershady, Jan- 
gren 2000; Bournaud et al. 2005), Hi distributions (e.g. 
Baldwin, Lynden-Bell, Sancisi 1980), integrated Hi line- 
profiles (Richtcr & Sancisi 1994), or kinematic features 
(e.g., Richter & Sancisi 1994; Swaters et al. 1999; Ko- 
rnreich et al. 2000, Chemin et al. 2006), especially in 
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LSB galaxies (Swaters et al. 2003; Simon et al. 2003). 
However, a direct measurement of disk ellipticity is chal- 
lenging because it is difficult to disentangle the intrin- 
sic shape of galaxy disks from both projection effects 
and spiral structure (e.g. Rix & Zaritsky 1995; Schoen- 
makers et al. 1997, Schoenmakers 1999; Andersen et al. 
2001; Barnes & Sellwood 2003). If better estimates of 
ellipticity existed, tighter limits could be placed on halo 
triaxiality. 

This paper describes our disk structure survey. With 
this survey, we attempt to constrain the halo shapes of 
field galaxies by measuring disk lopsidedness and elliptic- 
ity using a combination of bi-dimcnsional spectroscopy 
and imaging, complementing other work such as the 
Fabry-Perot kinematic survey of Virgo cluster disk galax- 
ies (Chemin et al. 2006). Beyond the implications such 
asymmetry measurements have for cosmological simu- 
lations of galaxy halos, both quantities also should be 
sources of scatter in the TF scaling relation (Franx & 
de Zeeuw 1992; Zaritsky & Rix 1997; Andersen & Ber- 
shady 2003). Therefore, a measurement of ellipticity, and 
possibly lopsidedness, places limits on other astrophys- 
ical sources of TF scatter, including disk mass-to-light 
ratios. If limits are placed on the variation of disk mass- 
to-light ratios, a new constraint will be placed on mass 
decomposition models. 

In this paper we present the Ha and H I kinematic data 
for our survey. Ha echelle observations were made with 
the integral field unit (IFU) DensePak (Barden, Sawyer 
& Honeycutt 1998) on the WIYN 3.5m telescope 3 for 39 
nearly face-on disk galaxies. Hi observations for 23 of 
the sample come from the Nangay radio telescope. The 
optical data forms the primary kinematic data for this 
survey. The purpose of the Hi line-profiles is to estab- 
lish that the optical integral-field data amply samples the 
potential. While Hi and Ha have similar flux- weighted 
radial extents in many cases (Verheijen & Sancisi 2001), 
there will still exist different and varying spatial distri- 
butions and "turbulent" motions of ionized and neutral 
gas which could lead to significant differences between 
Hi and Ha integrated line- widths. However, our ex- 
pectations are that turbulent motions are small; the dif- 
ferences between ionized and neutral turbulent motions 
smaller still; and the differences in spatial distribution 
are mitigated by the typically flat, and steeply rising 
rotation curves of this sample. Verifying these expec- 
tations is important for our ultimate goal of measuring 
halo shapes, and also allows us to address the suitability 
of optical line- widths as dynamical tracers of total mass. 
The application has relevance to studies of both nearby 
and distant galaxies. The combination of optical and ra- 
dio data also is pertinent to understanding the origin of 
line-profile asymmetries. 

In future papers in this series we analyze the Ha and 
Hi line-profiles, the Ha velocity fields, and the photo- 
metric properties of the survey galaxies. These studies 
explore the photometric and kinematic asymmetries of 
the galaxy sample, and their correlations. The studies 
build a physical picture of why asymmetries arise; how 
they contribute to scatter in the Tully-Fisher relation; 

3 The WIYN Observatory is a joint facility of the University of 
Wisconsin-Madison, Indiana University, Yale University, and the 
National Optical Astronomy Observatories. 



and how well total dynamical mass can be measured 
in nearly face-on systems. Together, they allows us to 
study the ellipticity distribution of our sample in a ter- 
tiary analysis. 

The present paper is organized as follows. We provide 
the selection algorithm in §2. We present the DensePak 
and Nangay observations in §3. The basic processing of 
the data is detailed in §4. Hi line-widths are also pre- 
sented in §4. In §5 we describe how we combine integral- 
field observations to obtain Ha integrated line- widths, 
present measurements of turbulent motions, discuss in- 
tegrated line-profiles of individual galaxies, and contrast 
optical and radio line-profiles. We conclude with a sum- 
mary in §6. 

2. SELECTION ALGORITHM 

The selection of the sample for this survey is driven by 
our desire to measure disk ellipticity. A 10% disk elliptic- 
ity will create differences in observable quantities which 
are smaller than typical measurement errors for inclined 
galaxies (Barnes & Sellwood 2003; Andersen, Bershady, 
Sparke 2005), thereby necessitating that our sample in- 
clude only nearly face-on galaxies. To illustrate this, 
Figure 1 shows how measurement errors on observable 
quantities translate into intrinsic ellipticity errors as a 
function of inclination. We have assumed that intrinsic 
ellipticity is measured via a comparison of observed pho- 
tometric and kinematic position angles, photometric axis 
ratios, and kinematic inclinations. A description of the 
method is found in Andersen et al. (2001). Given the 
expected errors on these quantities and that the mean in- 
trinsic ellipticity of galaxies is expected to be below 10%, 
galaxies with inclinations < 30° are optimum choices for 
intrinsic ellipticity measurements. 

Two methods exist for making a selection of nearly 
face-on disks: Galaxies can be chosen with photometric 
axis ratios close to unity, or chosen with narrow H I line- 
widths (or other kinematic measures of rotation). The 
Tully-Fisher relationship can be inverted; line widths 
and absolute magnitudes can used to solve for inclina- 
tion (Rix & Zaritsky 1995). Because the scatter in the 
Tully-Fisher relation is small for normal spirals, these so- 
called "inverse Tully-Fisher" inclinations could be used 
to select galaxies with inclinations less than 30° . 

Choosing galaxies with existing Hi line-widths may 
bias a sample. Many galaxies are selected for Hi ob- 
servations on the basis of having photometric inclina- 
tions greater than 45° for inclusion in Tully-Fisher stud- 
ies. Those galaxies with narrow Hi line- widths may in- 
deed be nearly face-on, but in this case are necessarily 
highly elliptic. Hence surveys built upon the availability 
of H I line-widths as a selection criterion may preferen- 
tially include spiral galaxies with elliptical disks. We 
therefore choose galaxies for this study based on photo- 
metric axis ratios. 

However, disks with the greatest ellipticities will ap- 
pear nearly face-on only if the intrinsic minor axis lies 
near the inclination axis. This potential bias is a second- 
order effect which can be accounted for when modeling 
the distribution in ellipticities for the general popula- 
tion of spiral galaxies (Ryden 2004). We therefore chose 
galaxies with axis ratios near unity. 

The DensePak integral field unit on the WIYN tele- 
scope spans 30 arcsec x 45 arcsec, so targets were lim- 
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Fig. 1. — The 68% confidence limit error distribution in intrin- 
sic disk ellipticity, <x e , versus inclination. The lighter solid lines 
represent the contribution to the ellipticity error due to the mea- 
surement errors of 5°, 7°, 10° and 12° in the difference between 
kinematic and photometric position angles (cr^). The dashed lines 
are the contribution to the error on ellipticity due to 0.01, 0.03 and 
0.05 axis ratio measurement errors (<7(,/ a ). The dotted line is the 
ellipticity error due to kinematic inclination errors (<Tj). For each 
inclination angle, we used an empirical relationship between incli- 
nation and kinematic inclination to estimate the error (Andersen 
& Bershady 2003). The heavy solid line is the total ellipticity error 
as a function of inclination for the inclination error, = 10° and 
u b/a =0.03. Ellipticity errors were calculated assuming an elliptic- 
ity of 0.1 and that the intrinsic ellipticity major axis is aligned with 
the kinematic major axis. While the details of these curves vary 
slightly for different cllipticitics and phase angles between the el- 
lipticity and kinematic axes, it is always true that ellipticity errors 
increase with inclination. 



ited to diameters (-D25) between 45 and 75 arcseconds 
so they could be well sampled in a few pointings. Given 
this relatively small size, chosen galaxies will typically 
be some of the faintest and most distant galaxies in the 
catalogs of nearby galaxies. 

To maximize the number of galaxies available for study, 
we chose our sample from the LEDA database 4 which 
contains all galaxies in common catalogs (e.g. NGC, 
UGC, IC), plus additional galaxies only defined in the 
Principal Galaxy Catalog (PGC; Paturel et al. 1997). 
While the PGC is inclusive of many other galaxy cat- 
alogs, galaxies with diameters less than ~1 arcminute 
still fall below the completeness threshold (Figure 2). 
We required sample galaxies from LEDA to have: (1) 
axis ratio b/a + Ab/a > 0.87 (i - Ai < 30°) 5 , (2) ap- 
parent disk size 45" < D25 < 75", (3) galactic absorp- 
tion A B < 0.6 (Schlegel, Finkbeiner & Davis 1998), (4) 
declinations visible from WIYN, ($1950 > —10°, (5) ob- 
served B-band magnitudes, (6) observed recession veloc- 
ities, and (7) fc-types ranging from Sab to Sd (t-types 

4 http://leda.univ-lyonl.fr 

5 Here, Ab/a is the uncertainty in the measured axial ratio b/a, 
and Ai is the corresponding uncertainty in inclination as tabulated 
in LEDA. 
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Fig. 2. — Distribution of diameters (D25; top panel), blue 
magnitudes (B; middle panel) and recession velocities (V; bot- 
tom panel) for galaxies in the PGC meeting the selection criteria 
in §2.2 (filled circles) and galaxies observed as part of this sur- 
vey (open triangles). Specifically, these selection criteria include: 
(1) 1.5 < t-type < 8.5, (2) A B < 0.6, (3) <5i 95 o > -10°, (4) 
b/a.25 + Afo/a25 > 0.87. Completeness criteria for galaxies in the 
PGC are met for D25 > 0.8' and B < 15.5. Despite the incomplete- 
ness of the parent sample, we endeavored to choose a representative 
sample of disk galaxies which spanned a wide range in types and 
luminosities. 



between 1.5-8.5). We elected not to observe Sa galax- 
ies because these galaxies have less Ha emission than 
later types (Kennicutt, Edgar & Hodge 1989). Of the 
1300 targets which met these raw selection criteria, 753 
were observable during our longest WIYN DcnsePak runs 
(January and December 1999; h < RA < 13 h ). 

We obtained Second Generation Digitized Sky Survey 
(DSS2) images to further refine our sample selection. 
After visual inspection, we eliminated from our sample 
galaxies with bars, rings, interacting companions, or fore- 
ground stars that contaminated the inner isophotes. This 
step removed 60% of the sample from consideration. We 
measured scale lengths, position angles and axis ratios 
from the DSS2 images using the IRAF ellipse routine, 
and we further required (1) a relatively constant posi- 
tion angle in the outer portions of the disk, therein re- 
moving most warped galaxies, and (2) our measurements 
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Fig. 3. — The distribution in type (top panel) and luminosity 
(bottom panel) for galaxies in the PGC meeting the selection cri- 
teria in §2.2 (filled circles) and galaxies observed as part of this 
study (filled triangles). The selection criteria applied to sample 
of galaxies in the top panel include, (1) measured B-magnitudcs 
and (2) recession velocities, (3) A B < 0.6, (4) 5ig 50 > -10°, (5) 
b/a 25 + Ab/a 25 > 0.87, and (6) 0.'6 < D 25 < 2'. The bottom 
panel include all galaxies in the top panel which have 1.5 < t-type 
< 8.5. We broadened the sample by both type and luminosity 
upon completion of our pilot study (open squares; Andersen et al. 
2001). 



oib/a > 0.87 on DSS2 images at three scale lengths. Be- 
tween 0' 1 < RA < 13 h only 70 targets remained. As a 
test of the completeness of our selection criteria, we ex- 
amined the DSS2 images of every PGC galaxy in a 1 hour 
slice of right ascension and ascertained that our selection 
criteria successfully yielded every galaxy with a photo- 
metric inclination less than 30° in this slice of right as- 
cension. Through the course of our observations, we ob- 
served 32 of these 70 targets, plus 7 other galaxies chosen 
from the same parent sample of 1300 targets, but which 
lie between I3 h < RA < 2A h . The 39 observed galaxies 
are not representative of the parent sample, rather we 
attempted to span a large range of morphological types 
and luminosity (Table 1 and Figure 3). 

3. KINEMATIC DATA 
3.1. DensePak Ha Observations 
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Fig. 4. — Schematic of DensePak in the telescope focal plane. 
Fibers arc identified according to their position in the output slit. 
There are empty spaces in the slit between fibers 22-23, 26-27, 
33-34, 36-37, 41-42, 53-54 and 86-87. Solid black disks mark 
broken fibers in the fiber bundle. Spatial orientation is for a zero 
rotation offset. Fibers are numbered in order of live fibers along 
the reconfigured output slit. This numbering scheme differs slightly 
from the DensePak manual. 



We observed 39 sample galaxies with DensePak dur- 
ing 11 nights on 6 separate observing runs (Table 2). 
The DensePak IFU is a fiber optic-array mounted at the 
Nasmyth f/6.3 focus imaging port on the WIYN 3.5m 
telescope. The WIYN Observatory, located on Kitt Peak 
south of Tucson, Arizona, stands 6875 feet above sea level 
and is home to the second largest telescope on the moun- 
tain. Dedicated in October, 1994, the WIYN telescope 
employs an altitude-azimuth mount design. DensePak 
is an array of 91 fibers arranged in a seven by thirteen 
fiber rectangle sub-tending an area of 28" x 42". Fibers 
have an active core diameter of 2.81" (300 jitm). Cladding 
and buffer increase the total fiber diameter to 3.75" (400 
/Lim). According to the DensePak manual's astromctric 
diagram (Sawyer 1997, Figure 4) 6 the fiber spacing is 4.2 
arcsec (center-to-center). In addition to the 91 fibers ar- 
ranged in a rectangle, another 4 fibers are spaced around 
the rectangle roughly an arcminute from the center and 
are used to measure the "sky" flux. Of the original 91 
DensePak fibers, 86 fibers remain — 5 fibers are broken 
(Figure 4). 

This manual is currently available at 
www.noao.edu/wiyn/densepak.pdf. The astromctric diagram 
was subsequently published in Homeier & Gallagher (1999). Note 
there is some uncertainty in the adopted fiber spacing. The 
plate scale of 10 arcsec mm - 1 quote in the manual is incorrect, 
and should be 9.374 arcsec mm -1 for the bare RC focus of 
the WIYN port where DensePak was used (see, for example: 
www.noao.edu/wiyn/wiynfacts.html). This reconciles the differ- 
ent fiber diameters quoted here versus the manual and Homeier 
& Gallagher. However, it is unclear if the manual's astromctric 
diagram is based on this erroneous plate scale. For the purpose of 
this work we adopt the implied fiber spacing given by the numbers 
in the astromctric diagram. 
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TABLE 1 

Properties of survey galaxies tabulated from the PGC a 
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1.2 


0.85 


15.3 


0.08 


3152 


-18.0 


38908 


UGC 7208 


12 09 56.8 


+39 23 20 


Sbc 


4.4 


1.0 


0.91 


15.0 


0.13 


7078 


-20.1 


39728 


NGC 4275 


12 17 21.8 


+27 53 54 


Sb 


2.8 


0.8 


0.90 


14.0 


0.09 


2317 


-18.6 


46767 


NGC 5123 


13 20 58.7 


+43 20 50 


Sc 


5.9 


1.1 


0.89 


13.7 


0.06 


8323 


-21.6 


49906 


NGC 5405 


13 58 40.3 


+07 56 35 


Sc 


4.8 


0.8 


0.94 


14.9 


0.17 


6922 


-20.0 


55750 


IC 132 


15 37 53.7 


+20 50 28 


Sc 


5.4 


1.0 


0.92 


14.5 


0.27 


4525 


-19.8 


56010 


MCG 6-35-5 


15 44 37.1 


+33 22 38 


SBc 


6.2 


0.8 


0.96 


15.3 


0.14 


4468 


-18.8 


57931 


UGC 10357 


16 20 18.3 


+40 33 50 


Sbc 


3.5 


0.8 


0.88 


15.1 


0.05 


9280 


-20.5 


58410 


UGC 10436 


16 29 24.2 


+41 15 44 


SBc 


5.3 


1.0 


0.97 


14.5 


0.05 


9059 


-21.0 


70962 


MRK 318 


23 15 05.9 


+13 43 38 


Sbc 


3.8 


0.7 


0.80 


14.2 


0.20 


4455 


-20.0 


71106 


NGC 7620 


23 17 37.3 
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Sc 
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1.0 


0.96 


13.8 


0.34 


9582 


-22.1 


72144 


UGC 12740 


23 39 22.1 


+23 32 15 


Sc 


6.0 


0.6 


0.95 


16.6 


0.27 


10521 


-19.4 


72453 


UGC 12784 


23 45 32.0 


+17 11 49 


SBbc 


3.6 


1.3 


1.00 


14.7 


0.30 


9952 


-21.2 



a Paturcl et al. (1997) 

^dc Vaucoulcurs type for the galaxy (do Vaucoulcurs et al. 1991). 

C Morphological type placed in a pseudo-numerical scale (Paturel et al. 1997). 

^Diameter measured from the isophote at a B-band surface brightness of /ig = 25 mag arcsec — 2. 
Axis ratio measured from the isophote at a B-band surface brightness of fj, g = 25 mag arcsec— 2. 

f 

A s is calculated using the galactic absorption law of Cardclli et al. (1989). 
^Absolute B — band magnitude adopts Hq = 75 km s — 1 Mpc -1 and heliocentric recession velocity corrected for peculiar motion associated with the local group and Virgo 
cluster. 



DensePak feeds the WIYN Bench Spectrograph, a 
fiber-fed spectrograph designed to provide low to 
medium resolution spectra. We used the Bench Spectro- 
graph camera (BSC) and 316 lines/mm echelle grating 
in order 8 to cover 6500A< AA < 6900A, with a disper- 
sion of 0.195 A/pix (8.6 km s _1 /pix) and an instrumental 
FWHM of 0.51 A (22.5 km s" 1 ). We required the highest 
resolution available with this instrument; fitting velocity 
field models to galaxies with observed rotation velocities 
of less than 100 km s" 1 required centroids with accura- 
cies of roughly 5 km s _1 . The BSC images the spectro- 
graph onto a T2KC thinned SITe 2048x2048 CCD with 
24 /xm-pixels. The spectra are aligned along the columns 
of the CCD. The chip has a read noise of 4.3 e~ and was 
used with the standard gain of 1.7 e~/ADU. The sys- 
tem throughput for this setup is roughly 4% estimated 
from comparing our mean spectral continuum fluxes to 
calibrated i?-band images (Bershady et al. 2005). 



Since galaxy rotation curves typically peak at roughly 
two scale lengths (Courteau & Rix 1999; Willick 1999) 
and Ha emission is only detected out to 3-4 scale lengths 
(Rubin, Waterman & Kenney 1999; Dale et al. 2001), 
our goal was to spatially cover galaxies out to 2.5-3.5 
scale lengths. Galaxies with DSS2 scale lengths less than 
six arcseconds only required one DensePak position to 
cover several scale lengths. However, to improve fits to 
the velocity field we required that two spatial resolution 
elements (fibers) be used to sample one disk scale length. 
For galaxies with scale lengths less than six arcseconds, 
two overlapping DensePak positions offset by a half-fiber 
diameter were used. For scale lengths between 6" and 
9", we used two DensePak positions offset to cover a 
45" x45" area, and for scale lengths between 9" and 12", 
we used three DensePak positions. The first two posi- 
tions were offset 27" from each other. For the third po- 
sition, DensePak was rotated, centered, and offset south 
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TABLE 2 
DensePak Observing Log 
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a The Run Identification identifies the Principal Investigator (PI) and observers 
for each run: (a) PI: Bcrshady (UW allocation), Observers: Bcrshady & Andersen; 
(b) PI: Bcrshady, Gallagher, Sparke, Wilcots (UW allocation), Observer: Ander- 
sen; (c) PI: Bcrshady, Gallagher, Sparke, Wilcots (UW allocation), Observers: 
Andersen &c G. Madsen; (d) PI: Bcrshady (UW allocation), Observer: Andersen; 
(e) PI: Andersen (NOAO allocation), Observer: Andersen; (f) PI: Bershady &c An- 
dersen (UW allocation), Observer: Andersen. 



of the first two positions. These three DensePak posi- 
tions covered almost one square arcminute. In order to 
observe the Ha line emission at three scale lengths and 
improve cosmic-ray rejection, we typically took two, ^30 
minute exposures at each position. 

To center DensePak on our target we used a slit- 
viewing camera built into the fiber-mounting module 
known as WIFOE. Use of this camera involves insert- 
ing a pellicle for simultaneous viewing of the source and 
the rear-illuminated IFU. In this way, we aligned fiber 
45 on a nearby guide star, and then zeroed the tele- 
scope offsets. After moving to our target, we could see 
the target galaxy in the slit-viewing camera only if the 
galaxy had a high central surface brightness and the sky 
was dark and clear. If we saw the galaxy, we confirmed 



it was centered on fiber 45 and performed the first off- 
set. If the galaxy was not visible, we performed a blind 
offset. After guiding was established, we began obser- 
vations. Typically, after two thirty minute exposures, 
guiding was turned off and the second set of telescope 
offsets were applied. After guiding was re-established, 
we took our next two exposures. If a third position was 
required, we re-acquired the nearby guide star and ro- 
tated DensePak. Since DensePak was not on the rotator 
axis, the star shifted a few arcseconds as DensePak was 
rotated 90°. After re-centering the guide star on fiber 45, 
we moved DensePak back to the source, applied the final 
telescope offset and took our final two exposures after re- 
establishing guiding. We found that rotating DensePak 
added up to ten minutes of overhead to our observations; 
we therefore preferentially chose galaxies which required 
only two DensePak footprints to sample the disk out to 
2.5 scale lengths. 

We encountered problems with accurately positioning 
DensePak on-sky due to problems with the pellicle and 
non-guided offsets. For some of our observations, the 
pellicle did not properly align images in the slit-viewing 
camera, so our spectra were not centered as expected. 
However, when fitting velocity field models to multiple 
pointings, model centers were easily determined from the 
data. Hence, the unaligned pellicle did not diminish the 
quality of the data. However, inaccurate telescope offsets 
did complicate data handling. We performed tests at the 
telescope on the accuracy of our offsets by removing each 
offset in turn and re-centering on the guide star. We dis- 
covered that the guide star moved as much as a fiber 
diameter (3 arcseconds) after removing the telescope off- 
sets. Therefore, we are forced to let telescope offsets be 
a free parameter in our model fits, thereby introducing 
extra fit parameters which could have been avoided if 
reliable offsets were possible. More tests are needed to 
determine an accurate method for applying telescope off- 
sets using DensePak on WIYN. 

At the beginning and end of each night, dome flats 
and bias frames were taken. We also observed Thorium 
Argon (ThAr) line lamps at the beginning of the night, 
before and after the dewar was filled, and again at the 
end of the night. The BSC dewar has a short hold time 
of only 12 hours, and typically was refilled during the 
course of the night. We observed sub-pixel shifts of ~ 0.1 
pixels in the centroids of the ThAr emission-lines taken 
at different times. The largest shifts are seen between 
frames taken before and after dewar fills. When reduc- 
ing spectra, we calibrated wavelengths using the ThAr 
frame taken closest to a given observation, within a pe- 
riod between filling of the dewar. 

The quality of data from DensePak varied across the 
slit formed by the 90 fibers. In particular, comparisons 
of dome flat fluxes showed the first eight fibers in the top 
of the slit had less than half the flux of most fibers (see 
Figure 3 of Bershady et al. 2005). After field flattening, 
this low throughput was translated into greater values 
for the continuum noise in these fibers. We believe this 
poor performance may be due to: (1) fibers 1-8 undergo 
the tightest bend near the slit end of DensePak of all the 
fibers. Bershady et al. (2004) present evidence that this 
tight bend radius increases the focal ratio degradation of 
fibers for a similar integral field unit, SparsePak, thereby 
leading to greater losses in the spectrograph; (2) fibers 1- 
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8 are at the top end of the DensePak slit. If the slit is not 
centered within the baffles at the spectrograph entrance, 
these fibers could suffer from further vignetting; (3) al- 
ternatively, these fibers may be otherwise over-stressed, 
e.g. because these fibers are on an outer face of the array 
(see Bershady et al. 2004). Data from fiber 1 were so 
poor that we did not include them when, for example, 
we stacked spectra to make synthetic line profiles (§5.2). 

3.2. Nangay Hi Observations 

The Nangay telescope is a meridian transit-type in- 
strument of the Kraus/Ohio State design, consisting of 
a fixed spherical mirror, 300 m long and 35 m high, a 
tiltable flat mirror (200x40 m), and a focal carriage mov- 
ing along a 90 m long curved rail track, which allows the 
tracking of a source on the celestial equator for about I 
hour. Located in the center of France, it can reach de- 
clinations as low as -39°. It has an effective collecting 
area of roughly 7000 m 2 (equivalent to a 94-m diameter 
parabolic dish). Due to the elongated geometry of the 
mirrors, at 2I-cm wavelength it has a half-power beam 
width (HPBW) of 3'.6 E-Wx 22' N-S for declinations be- 
low 30°; at higher declinations the N-S HPBW increases 
(see plot in Matthews, van Driel & Monnier-Ragaigne 
2001). Typical system temperatures were ~40 K for our 
project. For a technical description of the Nangay deci- 
metric radio telescope and the general methods for data 
handling and reduction see Theureau et al. (1998a) and 
references therein. 

Observations at Nangay of 23 sample galaxies were 
made in the periods March to August 1999 and January 
to June 2001 using a total of about 110 hours of tele- 
scope time. We obtained our observations in total power 
(position-switching) mode using consecutive pairs of two- 
minute on- and two- minute off-source integrations. Off- 
source integrations were taken at approximately 20' E of 
the target position. The autocorrelator was divided into 
two pairs of cross-polarized receiver banks, each with 512 
channels and a 6.4 MHz bandpass. This yielded a chan- 
nel spacing of 2.64 km s _1 , for an effective velocity res- 
olution of ^3.3 km s _1 at 21-cm, which was smoothed 
to a channel separation of 7.91 and a velocity resolution 
of 9.50 km s _1 during the data reduction, in order to 
search for faint features. The center frequencies of the 
two banks were tuned to the known redshifted Hi fre- 
quency of the target. 

4. BASIC PROCESSING 

4.1. DensePak Reductions 

Data were overscan and bias-corrected and trimmed 
using the NOAO IRAF package cedproc. Since two 
or more frames were taken at each position, we cal- 
culated pair-wise differences and statistical thresholds, 
pixel- by-pixel, of the expected variance (a 2 ) due to read- 
noise and shot-noise. Pixels which were more than 
5cr brighter than the corresponding pixels on the other 
frames were flagged as cosmic rays and masked from the 
final, mean-combined, two-dimensional spectral image. 
This is preferable to removing cosmic rays after spectral 
extraction which averages cosmic rays falling in pixels 
that are part of extracted apertures with good data from 
other pixels in same spectral channel. Following cosmic- 
ray cleaning, basic spectral extraction, flattening, wave- 
length calibration and sky subtraction were done using 



IRAF dohydra. Extracted one-dimensional spectra were 
field-flattened with dome flats and wavelength-calibrated 
using ThAr emission spectra. Finally, the four sky spec- 
tra were averaged and the mean sky spectrum was sub- 
tracted from the other 86 spectra. 

Once spectra were processed with IRAF dohydra, 
we identified extracted Ha emission-line characteristics. 
Our algorithm measured Gaussian fluxes, widths, cen- 
ters and centroid errors for lines in a given spectral in- 
terval. Most emission-lines in our data were symmet- 
ric and were well-fit by a single Gaussian. Some Ha 
profiles sampled by fibers within a fiber radius of the 
galaxy center were skewed or even bimodal as would be 
expected when the fiber diameter is larger than the lo- 
cal dynamical scale. Since the fraction of such lines was 
small, we flagged these lines during visual inspection af- 
ter running our automated line fitting algorithm and do 
not cite the fit results later. While Beauvais & Bothun 
(1999) found the best fits to a variety of different sim- 
ulated line-profiles were obtained with Vogt profiles (a 
combination of Gaussian and Lorentzian profiles), they 
noted that Gaussian fits yielded comparable results for 
signal-to-noise (S/N) ratios greater than 20. A similar 
conclusion was reached by Courteau (1997) who mea- 
sured intensity-weighted centroids. Our typical Ha lines 
had high S/N, and a visual inspection of the fits made 
it clear that Gaussian profiles were sufficient to describe 
most emission lines. 

Since the S/N of each channel in the spectra is im- 
portant for accurate profile fitting, we iteratively estab- 
lished measurement errors for each channel. We begin 
by calculating an iterative clipped mean and standard 
deviation within a 30A window containing the emission 
line. We used a very tight clipping criterion: All channels 
with counts greater than the standard deviation were re- 
jected, repeating until the number of channels did not 
decrease further. This procedure removed almost all ef- 
fects of emission-lines, imperfectly subtracted sky lines, 
and residual cosmic rays from the continuum mean and 
standard deviation. We normalized the standard devia- 
tion of the iteratively clipped data by the standard de- 
viation of a Gaussian distribution which also had been 
truncated at la. With the continuum standard devia- 
tion, (Tcont., established, we assumed remaining errors on 
channel counts were based on Poisson counting errors. 
Therefore the error on each channel, <7,, was calculated 
as follows: 

u i = {"font. + (Counts, - Counts cont .)} 1/2 (1) 

where all counts were calculated in e~ /pixel, Counts^ 
was the counts for channel i, and Counts cont . was the 
iterative mean calculated in the 200A window. 

Once electron counts and the initial estimates of the 
standard deviations had been calculated for every chan- 
nel, we used the Marquardt-Levinson algorithm to fit a 
Gaussian profile plus a linearly varying background to 
the line within the selected region. We then subtracted 
this best fit spectrum within the window and recalculated 
the standard deviation following the same procedure de- 
scribed above. Then we refit the data and extract the 
Gaussian amplitude, center (Ac) and width (ac) from 
the fit. The area beneath the curve yields the flux (/<?)■ 
The error on the centroid (a\ c ) was obtained from the 
covariance matrix. 
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Fig. 5. — Ho emission-linc-profilcs for 6 fibers selected from the first pointing of PGC 14564. Fiber 30 was a sky fiber and illustrates 
the typical background fluctuations. The solution for a Gaussian Ha emission-line in Fiber 24 (upper right panel) did not meet the S/N 
criterion: £ = /G/°"cont./(2o" c v) 1 / 2 > 5. For the emission-line in aperture 10 (left, middle panel), f = 8 . The Ha emission-line for aperture 
25 (middle right panel) was typical for most observed galaxies. Aperture 27 (bottom left) shows very strong Ha emission for our survey. 
The Ha profile in aperture 47 (lower right panel) was highly skewed, showing signs of bimodality. 
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We found that features which had Gaussian fits that 
did not reach the S/N threshold /G/o'cont./( 2(7 G) 1/ ' 2 > 5 
were also not identified as emission-lines by visual in- 
spection. We also required the width of the line must be 
equal to or greater than the resolution element. Finally, 
as already noted, we discarded lines that, based on vi- 
sual inspection, were skewed, required multiple Gaussian 
profiles, or formally met the S/N and width requirements 
but none-the-less did not look plausible. Figure 5 shows 
a sample of the output with Gaussian fits chosen for their 
range of properties: an example of sky continuum, a line 



which did not meet the S/N threshold, a weak line which 
did have sufficient S/N, typical, strong and skewed lines. 

Table 3 contains a sample of the DensePak line-fits for 
the first pointing of PGC 14564. The line-fits tables for 
all 39 galaxies are available in electronic form. Tabulated 
offsets refer to the fiber positions West and North of the 
kinematic center based on our model fits to the velocity 
fields described in the relevant forthcoming paper in this 
series, and the astrometric information regarding relative 
fiber positions described in §3.1. 



TABLE 3 

Gaussian Fits to Ha emission-lines: 1st pointing of PGC 14564 
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9 


24 


1 


-18.43 


-12.67 


6638.57 


0.17 


331 


0.52 


9 


25 


1 


-14.28 


-12.67 


6637.82 


0.02 


2411 


0.51 


8 


26 


1 


-10.13 


-12.67 


6637.75 


0.01 


9654 


0.50 


10 


27 


1 


-5.98 


-12.67 


6637.72 


0.01 


16032 


0.49 


12 


28 


1 


-1.83 


-12.67 


6637.64 


0.01 


12010 


0.51 


12 


29 


1 


2.32 


-12.67 


6637.84 


0.01 


9437 


0.50 


14 


30 


1 


-56.3 


42.5 










8 


31 


1 


-20.50 


-9.07 


6638.73 


0.14 


527 


0.64 


9 


32 


1 


-16.35 


-9.07 


6638.54 


0.04 


960 


0.50 


8 


33 


1 


-12.20 


-9.07 


060 / .95 


0.02 


4054 


0.58 


9 


34 


1 


-3.90 


-9.07 


6637.86 


0.01 


17991 


0.50 


15 


35 


1 


0.25 


-9.07 


6637.87 


0.01 


12811 


0.51 


17 


36 


1 


4.40 


-9.07 


6638.19 


0.01 


12383 


0.63 


13 


37 


1 


-18.43 


-5.48 


6638.84 


0.03 


1256 


0.37 


8 


38 


1 


-14.28 


-5.48 


6638.81 


0.02 


3654 


0.50 


8 


39 


1 


-10.13 


-5.48 


6638.46 


0.01 


8681 


0.62 


10 


40 


1 


-5.98 


-5.48 


6638.20 


0.01 


17109 


0.60 


16 


41 


1 


-1.83 


-5.48 


6638.19 


0.01 


17713 


0.64 


21 


42 


1 


-20.50 


-1.89 


6638.88 


0.06 


598 


0.41 


8 


43 


1 


-16.35 


-1.89 


6638.87 


0.01 


4132 


0.46 


9 


44 


1 


-12.20 


-1.89 


6638.89 


0.01 


8307 


0.51 


9 


45 


1 


-8.05 


-1.89 


6638.80 


0.01 


13259 


0.67 


11 


46 


1 


-3.90 


-1.89 


6638.73 


0.01 


17517 


0.76 


22 


47 


1 


0.25 


-1.89 


6639.10 


0.02 


21608 


1.04 


36 


48 


1 


4.40 


-1.89 


6639.50 


0.01 


19891 


0.97 


28 


49 


1 


-22.6 


1.7 










8 


50 


1 


-18.43 


1.71 


6639.04 


0.04 


667 


0.37 


7 


51 


1 


-14.28 


1.71 


6639.17 


0.01 


5098 


0.61 


8 


52 


1 


-10.13 


1.71 


6639.34 


0.01 


10843 


0.76 


10 


53 


1 


-5.98 


1.71 


6639.78 


0.01 


17740 


0.96 


17 


54 


1 


2.32 


1.71 


6640.30 


0.01 


18074 


0.92 


32 


55 


1 


-20.5 


5.3 










8 


56 


1 


-16.35 


5.30 


6639.72 


0.03 


2233 


0.59 


9 


57 


1 


-12.20 


5.30 


6639.93 


0.01 


11948 


0.60 


10 


58 


1 


-8.05 


5.30 


6640.27 


0.01 


19399 


0.73 


12 


59 


1 


-3.90 


5.30 


6640.67 


0.01 


24179 


0.68 


17 


60 


1 


0.25 


5.30 


6640.92 


0.01 


17069 


0.67 


20 
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TABLE 3 — Continued 



Fiber Pointing West North Xq AAp fa °"G o"cont. 

(arcscc) (arcscc) (A) (A) (e _1 /hr) (A) (e _1 /hr/res) 



61 


I 


4.40 


5.30 


6641 05 


0.01 


8105 


0.57 


16 


62 


I 


-57.3 


-48.5 










g 


63 


1 


-22.6 


8.9 










9 


64 


1 


-18.4 


8.9 










9 


65 


1 


-14.28 


8.90 


6640 03 


0.01 


7204 


0.55 


9 


66 


1 


-10.13 


8.90 


6640 19 


0.01 


14114 


0.59 


10 


67 


1 


-5.98 


8.90 


6640 78 


0.01 


20675 


0.68 


12 


68 


1 


-1.83 


8.90 


6641.15 


0.01 


20157 


0.63 


17 


69 


1 


2.32 


8.90 


6641.42 


0.01 


10514 


0.55 


16 


70 


1 


-20.5 


12.5 










10 


71 


1 


-16.4 


12.5 










9 


72 


1 


-12.20 


12.49 


6640 40 


0.02 


3521 


0.58 


10 


73 


1 


-8.05 


12.49 


6640.74 


0.01 


7907 


0.61 


12 


74 


1 


-3.90 


12.49 


6641.18 


0.01 


11172 


0.61 


11 


75 


I 


0.25 


12.49 


6641 ^5 

Uut:1 . <JtJ 


0.01 


8510 


0.53 


14 


76 


I 


4.40 


12.49 


6641 66 


0.01 


9271 


0.51 


13 


77 


I 


-22.6 


16.1 










9 


78 


1 


-18.4 


16.1 










10 


79 


1 


-14.3 


16.1 










10 


80 


1 


-10.13 


16.08 


6640.89 


0.02 


2076 


0.49 


9 


81 


1 


-5.98 


16.08 


6641.10 


0.01 


5295 


0.53 


11 


82 


1 


-1.83 


16.08 


6641.40 


0.02 


5337 


0.57 


14 


83 


1 


2.32 


16.08 


6641.69 


0.01 


7762 


0.54 


11 


84 


1 


33.7 


-54.5 










8 


85 


1 


-20.5 


19.7 










13 


86 


1 


-16.4 


19.7 










10 


87 


1 


-8.05 


19.68 


6641.19 


0.02 


2554 


0.50 


11 


88 


1 


-3.90 


19.68 


6641.26 


0.01 


5194 


0.49 


13 


89 


1 


0.25 


19.68 


6641.40 


0.02 


3447 


0.49 


12 


90 


1 


4.40 


19.68 


6641.64 


0.03 


2670 


0.53 


12 



We followed this procedure for the sample of 39 galax- 
ies, or 8280 individual spectra. 6000 (72.5%) of the spec- 
tra showed clear emission lines that were detected and 
well fit by the algorithm based on visual inspection. 1875 
(22.6%) of the spectra had no detectable emission lines. 
368 of these were sky spectra. Visual inspection nagged 
46 (0.6%) spectra with potential emission lines that were 
undetected by the algorithm. All of these had low S/N 
and indeed simply may not have met the S/N criteria. 
98 (1.2%) of the spectra were asymmetric, showing ei- 
ther a high degree of skew or requiring multiple compo- 
nents. A single Gaussian profile is not a good fit to these 
profiles; widths and centroids are quite possibly in error. 
They were excluded from Table 3. As stated above, these 
fibers are located very close to the center. In most cases, 
they are within a fiber radius of the center. 48 (0.6%) of 
the spectra had fits but were ruled out by visual inspec- 
tion. Most of these false positives are weak cosmic rays 
or a result of some residual curvature in the spectra that 
was not removed through flat-fielding. These are not in- 
cluded in Table 3. The remaining 213 (2.6%) of spectra 
had questionable fits. These algorithm detections are at 
a low S/N, and may not have been detected as lines by 
visual inspection. Some of these questionable fits may 
be cosmic rays. Most are just emission lines which just 
barely meet the detection threshold. 

4.2. Hi Reductions 

We reduced our H I spectra using the standard Nangay 
spectral line reduction packages available at the Nancay 
site. With this software we subtracted baselines (gener- 
ally third order polynomials), averaged the two receiver 
polarizations, and applied a declination-dependent con- 
version factor to convert from units of T svs to flux den- 



sity in mJy. The T sys -to-mJy conversion factor is deter- 
mined via a standard calibration relation established by 
the Nancay staff through regular monitoring of strong 
continuum sources. This procedure yields a calibration 
accuracy of ~15%. In addition, we applied a flux scaling 
factor of 1.26 to our spectra based on statistical compar- 
isons (see Matthews et al. 1998, Matthews et al. 2001) of 
Nancay survey data of samples of late-type spirals with 
past observations of these galaxies made at Nancay and 
elsewhere. 

From the 23 galaxies observed with Nancay from the 
two campaigns (March- August, 1999 and January- June 
2001), 18 had sufficient signal to measure line-widths and 
integrated Hi fluxes (Table 4). The linewidths listed 
in Table 4 have not been corrected for redshift stretch. 
Restframe linewidths should be reduced by a factor (1 + 

4.3. Hi Values from the Literature 

Hi observations from the literature were drawn from 
the LEDA extragalactic database (Paturel 2003) to com- 
plement Nancay observations. Corrections were made to 
the raw H I data to account for velocity resolution and in- 
homogeneity of the various original references to produce 
velocity widths measured on a common system using the 
optical convention (v = cAX/X) at 20% and 50% of the 
peak flux. These corrected W20 and W50 line widths are 
tabulated for 28 sample galaxies in Table 5. Hi obser- 
vations either from the Nancay observations described 
above or from the LEDA database exists for 38 of the 
39 sample galaxies. There are 10 galaxies with W50 line 
widths from both Nancay observations and the litera- 
ture, and the measurements are in good agreement. The 
mean difference in Nangay versus LEDA W50 line widths 
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TABLE 4 

H i Parameters from Nanqay Observations 



PGC 


rms a 


S/N b 


Ihi c 


w 50 d 


W 20 e 


V H i { 


* HI 


log M HI h 




( m Jv) 




(Jv km s — M 
v u j 11111 " ) 


(km s" 1 ) 


(km s" 1 ) 


(km s" 1 ) 


(km s" 1 ) 


floe Mr>) 


03512 


5.82 










5436 






08941 


5.31 


6.69 


4.72 


162 




9447 


5.6 


10.16 


14564 


2.96 


26.33 


10.06 


176 


202 


3483 


2.2 


9.71 


15531 


5.60 


5.63 


4.23 


157 


197 


5552 


12.9 


9.74 


16274 


3.74 










8883 






19767 


4.21 


10.34 


5.35 


141 


164 


5136 


5.4 


9.77 


23598 


2.80 


10.41 


2.57 


119 


139 


7478 


4.9 


9.78 


23913 


3.52 


34.82 


4.39 


33 


57 


6330 


1.6 


9.87 


24788 


2.86 


11.56 


4.63 


161 


194 


7545 


5.8 


10.04 


26140 


1.79 










8760 






28310 


2.7 


7.55 


2.65 


141 


153 


5905 


5.1 


9.59 


31159 


4.9 


3.26 


2.32 


190 


199 


5781 


10.7 


9.51 


32091 


3.03 










2511 






32638 


2.80 


5.45 


1.83 


172 


229 


6527 


15.8 


9.51 


33465 


4.41 










5935 






36925 


7.67 










6826 






38908 


3.87 


13.84 


7.61 


195 


221 


7084 


4.2 


10.20 


39728 


1.97 


8.81 


2.08 


145 


196 


2316 


9.1 


8.67 


46767 


2.89 


8.47 


4.06 


276 


335 


8065 


10.4 


10.07 


56010 


4.80 


4.78 


4.80 


63 


160 


4469 


23.2 


9.60 


57931 


6.09 










9280 






58410 


5.21 


5.09 


4.44 


256 


326 


9055 


18.9 


10.18 


72453 


3.24 


11.26 


4.38 


160 


194 


9932 


6.0 


10.26 



a Root mean squared noise level of the spectrum. 

^Ratio of peak line flux to the rms noise level. During reduction, data was smoothed to an effective velocity resolution of 9.5 km/s. 
c Intcgrated line flux. 

^Velocity width at 50% the peak line flux. According to Fouque et al. (1990), the uncertainty in W50 is 2av HI - 
c Vclocity width at 20% the peak line flux. According to Fouque et al. (1990), the uncertainty in W20 is 3(Tv HI ■ 
^Heliocentric central radial velocity of a line profile in the optical convention. 

g Uncertainty in Vhi, is a v HI — 4i? 1 / 2 P^ 2 (S/./V) _ 1 where R is the instrumental resolution (15.6 km s _1 ), Pw — (W20 — W50)/2 (Fouque et al. 
1990). 

h Total Hi mass, M HI — 2.356 x 10 4 D 2 I H /where D is the distance calculated assuming H — 75 km s" 1 Mpc -1 . 



is 1.7 km/s; the greatest difference in these W50 values 
(15.3±12.6 km/s) is consistent with the errors. Three 
of the six galaxies with both Nangay and literature W20 
line widths show discrepancies greater than 20 km/s. In 
these cases, the relatively low signal to noise ratios of 
the Nangay observations for PGC 39728, PGC 46767, 
and PGC 58410 have led to an over-estimate of W20 (see 
§5.4.1). 

5. Ha LINE- WIDTHS 

Creating Ha linewidth profiles from DensePak obser- 
vations is a computationally straight-forward process, 
as described below, but it is necessary to ensure these 
linewidths are physically meaningful. The filling fac- 
tor of DensePak is only 45% (factoring in the 4 dead 
fibers plus the highly vignetted fiber 1). In addition to 
the low filling factor, Ha may exist at radii beyond the 
three scale lengths typically covered by the DensePak ob- 
servations. Further, since the sample is nearly face-on, 
the "turbulent" component of the spatially-integrated 
line-width is critical to understand. Therefore, we first 
characterize the spatially-resolved Ha line- widths, as ob- 
served with DensePak. After describing the construction 
of the DensePak Ha linewidth profiles, we compare the 
DensePak Ha linewidth profile to one generated from 
narrow-band Ha imaging. We also comment briefly on 
the shape and width of the Hi and Ha profiles, with a 
more detailed analysis presented in the relevant forth- 
coming paper in this series. 



5.1. Spatially-Resolved DensePak Line-widths 

The projected Ha velocity widths from individual 
fibers vary little with radius, when considering the sam- 
ple taken as a whole. Variations do exist, however, be- 
tween and within galaxies; these variations are explored 
in a future paper. The sample-mean near-constancy with 
galactic radius is similar to what has been found for indi- 
vidual, nearly face-on disk galaxies in neutral gas (NGC 
628, NGC 1058 and NGC 3938; van der Kruit & Shostak 
1982, 1984; Dickey et al. 1990), molecular gas (NGC 628 
and NGC 3938; Combes & Becquaert 1997), and ionized 
gas (NGC 3938 and NGC 5668; Jimenez- Vicente et al. 
1999; Jimenez- Vicente & Battaner 2000). The reasons 
for this gross similarity of the near-constancy with ra- 
dius of different phases of the ISM is outside the scope of 
this paper, but we do note that some radial dependence 
in the neutral gas is observed in some galaxies (Shostak & 
van der Kruit 1984; Kamphuis & Sancisi 1993), and vari- 
ations are also seen in the ionized gas in some cases in our 
sample, albeit over different radial ranges. Our distribu- 
tion median, out of 6000 fiber samples with emission lines 
well characterized by a single Gaussian across 39 galax- 
ies in our sample is 17.9±2;g±|;Y km s _1 , where the two 
sets of "errors" include 50% and 80% of the distribution. 
These values are corrected for instrumental broadening 
only. They are consistent in value and distribution for 
what has been observed in two other, nearby spirals NGC 
3938 and NGC 5668 (Jimenez- Vicente et al. 1999, 2000), 
and are roughly two to three times the dispersion mea- 
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TABLE 5 

Hi Line-Widths from Paturel et al. (2003) 



PGC 


w 20 


°"w 20 


w 50 


°"W50 , 




(km s- 1 ) 


(km s- 1 ) 


(km s- 1 ) 


(km s- 1 ) 


03512 


204.9 


10.7 


182 


8.1 


05345 


52.49 


13.32 


27.05 


9.7 


05673 


147.75 


12.85 


134.05 


9.7 


06855 


136.5 


10.7 


115.7 


8.1 


07826 


96.08 


8.32 


85.4 


8.24 


08941 


157.22 


12.85 


147.32 


6.86 


14564 


198.93 


7.39 


179.64 


4.62 


19767 


147.9 


15.86 


137.14 


11.94 


20938 






169.61 


8.5 


23333 


139.12 


13.28 


134.42 


7.1 


23598 






121.37 


10.49 


23913 






47.37 


12.51 


24788 






147.76 


17.61 


26517 


126.5 


10.7 


125.7 


8.1 


27792 






129.87 


8.47 


28401 


59.5 


10.7 


37.7 


8.1 


33465 


285.9 


15.86 


239.14 


11.94 


3fiQ95 


149.9 


15.86 


138.14 


11.94 


38268 


183.33 


11.97 


162.92 


10.02 


39728 


163.28 


7.08 


148.96 


7.46 


46767 


292.9 


15.86 


253.14 


11.94 


55750 


108.05 


6.41 


92.61 


16.26 


57931 


176.41 


12.54 


170.19 


11.78 


58410 


248.9 


15.86 


231.32 


8.68 


70962 


242.35 


6.03 


209.09 


9.73 


71106 


266.78 


8.46 


236.86 


7.34 


72144 






125.31 


6.97 


72453 


179.2 


13.28 


157.3 


10.02 




10 20 30 40 



a„ a (km/s) 

Fig. 6. — Distribution of Ha emission-line velocity widths (a), as 
observed for the entire sample of 39 galaxies, corrected for instru- 
mental broadening. The instrumental resolution limit and integral 
distribution percentiles are indicated. The instrumental resolution 
is well below the peak of the observed distribution. 

spectra from fiber 1, as they suffered from extremely low 
signal to noise. 



sured in either neutral (van der Kruit & Shostak 1982, 
1984) or molecular gas (Combes & Becquaert 1987) for 
relatively small samples of galaxies. The tightness of our 
sample distribution is shown in Figure 6. The largest 
values which skew the distribution slightly are seen pri- 
marily at the center of galaxies with steep inner rotation 
curves. While beam smearing appears to affect these 
values, the effect is small, correctable, and does not ad- 
versely influence our results here. 

5.2. Spatially-Integrated DensePak Line-widths 

After taking the steps as described in §4.1 of extracting 
the spectra and fitting Gaussian line-profiles, we gener- 
ated "spatially integrated" line-widths - what we refer to 
here also as "synthetic line profiles." We coadded either 
the individual extracted spectra (channel by channel), or 
we coadded the individual fitted Gaussian line-profiles for 
those spectra where fits were possible. When combining 
these data or fitted profiles, we maintained a count of the 
number of independent fiber positions. An on-sky loca- 
tion was assigned to each fiber based on the DensePak 
geometry and telescope offsets determined via velocity- 
field fitting (presented in a future paper). The location 
of fibers from different pointings were checked for over- 
lap. If there was an overlap between fibers, the rele- 
vant fiber flux was scaled such that the combined spectra 
gave equivalent spatial sampling to non-overlapping ar- 
eas. Specifically, for each fiber we determined the number 
of fibers in a hexagonal area inscribed within a 6.3 arcsec 
radius circle centered on the fiber, and compared this to 
the expected number in non-overlapping regions (e.g., 7 
fibers for the middle of the array). We did not include 



5.3. Spatially-Integrated Line-widths from Narrow-band 

Imaging 

To test whether the Ha linewidth profiles derived from 
DensePak observations were valid representations of the 
matter and velocity distribution of the ionized gas, we 
constructed a linewidth profile from the narrow-band Ha 
image of one galaxy in our sample, PGC 38268, for com- 
parison. 

5.3.1. Observations 

PGC 38268 was observed in narrow-band filters using 
the SPIcam camera at the Apache Point Observatory 
(APO) 3.5m telescope on April 11, 2002. SPIcam is an 
imager equipped with a 2048 x 2048 pixel SITe back- 
illuminated CCD that has a 4.8' field of view. In the 
standard mode of operation (2 x 2-binning) SPIcam has 
an effective plate scale of 0.281"/pixel. Rest-frame Ha 
and a 6629A off-band filter come from a set owned by P. 
Hodge. SPIcam is designed to accommodate 3 x 3-inch 
filters, so the 2 x 2-inch narrow-band filters consequently 
reduce SPIcam's effective field of view (FOV) . 

The recession velocity of PGC 38268 shifted the Ha 
emission into the 6629A narrow-band filter. The Ha and 
6629A filters make an ideal pair of filters for narrow- 
band imaging because the filters' transmission curves are 
virtually the same in width and shape. Since the night 
was not photometric, no attempt was made to collect flux 
standards. Exposures in each band totaled 10 minutes. 

5.3.2. Reductions 

We used standard IRAF tasks to process images (over- 
scan correction, bias subtraction, and flat-fielding). Im- 
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ages were trimmed to include only the un- vignetted FOV. 
The standard IRAF cosmic ray cleaning task was aug- 
mented using information from neighboring pixels, thus 
enhancing the cleaning of extended cosmic ray events and 
detector flaws. 

Once individual data frames were fully processed, 
the sky background was calculated using an iteratively 
clipped mean and subtracted from all frames. On and 
off-band images were registered using 5 field stars. Since 
the ratio of fluxes in on versus off-band for these stars 
were in agreement with the predicted ratio based on the 
filter transmission curves, we used the ratio expected 
from the transmission curves to scale the images and sub- 
tract frames. The Ha narrow band image of PGC 38268 
with the overlay of DensePak fibers in which Ha was 
detected is shown in the first panel of Figure 7. While 
DensePak does not detect all the Ha flux because the the 
two clumps of Ha to the SW and the faint clumps to the 
N were not covered by the DensePak array, DensePak is 
able to detect very faint levels of Ha emission which are 
not visible in the image. 

5.3.3. A Comparison of Line Profiles 

We use this Ha narrow-band image of PGC 38268 to 
produce a linewidth profile. To do this, we assign a 
velocity to each position in the image. Following the 
fitting procedure outlined in Andersen et al. (2001), a 
best fit velocity field was determined from the DensePak 
data. The second panel of Figure 7 shows the observed 
velocity field from the DensePak data overlaid on the 
narrow-band image of PGC 38268. The best fit analytic 
model to this observed velocity field was then used to 
assign a velocity to each pixel of the narrow-band image. 
The resultant profile (dashed line of Figure 7) was then 
smoothed by a 20 km s _1 width Gaussian profile (solid 
line of Figure 7). This smoothing corresponds to the 
median measured Ha line- width in spectra of individual 
DensePak fibers. 

When the linewidth derived from the DensePak data 
if overlaid on this image (dotted line of Figure 7) , good 
agreement is observed between the smoothed profile and 
the DensePak profile. W20, the velocity width measured 
at 20% of the peak velocity, when measured from the 
profile derived from the Ha image is 152 km s -1 and 
when measured from the profile created directly from the 
DensePak spectra is 162 km s _1 . A difference of 5 km s _1 
in the width of the smoothing Gaussian used to smooth 
the line profile can account for this difference. The more 
striking discrepancy is between the peak flux levels in 
the two halves of the profile. The ratio of fluxes from the 
DensePak data is 0.6 while the ratio is 0.5 when mea- 
sured from the Ha image. This agrees qualitatively with 
the relative positions of Ha emission and the DensePak 
footprint: More Ha flux falls between fibers in the re- 
ceding (Northern) half of the velocity field. Finally, the 
difference in velocity between the two peaks is slightly 
greater in the data from the Ha imaging. This is con- 
sistent with the non-detection of Ha at large radii and 
velocities in the DensePak data. 

Overall, however, we conclude the DensePak-derived 
integrated line-width is likely an excellent approxima- 
tion. 

5.4. Synthetic Ha Line-Profiles 



We present these synthetic line-profiles for both the 
raw data and the summed, best-fit Gaussian profiles in 
Figure 8, along with an overlay of the DensePak obser- 
vations on the i?-band images. As can be seen by di- 
rect visual inspection, the line profile width and shape is 
very similar for the raw and Gaussian summations. The 
smoothing of the data with the Gaussion summation does 
not quantitatively alter low-order moments of the inte- 
grated line-profiles. For example, we find the variation of 
the skew, or lopsidedness, of the line profile between raw 
and Gaussion summations is less than the observational 
error, and are identical in the mean. 

A comparison between the H 1 and Ha profiles are also 
presented in Figure 8, where the data are available. In 
general there is very good agreement between the widths 
of the Hi and Ha profiles. Using the Gaussian fits to the 
individual Ha spectra, we constructed a width by finding 
the maximum difference in centroid velocities. We com- 
pared this width to the 31 Hi W20 measures and found a 
constant offset between the relations. The Hi line widths 
are 23±20 km/s broader than the Ha velocity "widths" 
as defined above, roughly comparable to the ~18 km/s 
spatially-resolved line-widths of the Ha emission lines. 
Only two of these galaxies, PGC 56010 and PGC 70962, 
showed significantly larger differences in velocity widths. 
As noted below, only half the disk of PGC 56010 is cov- 
ered with DensePak pointings. Similarly, we only cover 
the inner starbursting region of PGC 70962, but we sus- 
pect the width of the Hi profile may be dominated by 
the especially broad emission lines observed in the Ha 
spectra and may not reflect the kinematics of the galaxy. 
Further comparisons of H I and Ha line- widths are in a 
forthcoming paper. 

5.4.1. Notes to Individual Galaxies 

For galaxies with both DensePak Ha and Nangay 
Hi line-profiles shown in Figure 8, we normalized the 
peak fluxes of the Ha line-profiles generated from the 
sum of individual Gaussian fits to the Hi line profile 
peak flux for comparison. Here, we take note of the fea- 
tures observed in the line-widths for galaxies with both 
DensePak Ha and Nangay H 1 data: 

PGC 3512: No Hi flux from the source was detected. 

PGC 08941: Nangay observations were severely af- 
fected by radio interference, but no Hi flux was 
detected at the recession velocity of the source. 

PGC 14564: Hi and Ha line-profiles are very similar. 
The only notable difference is the slightly broader 
width of the Ha profile. 

PGC 15531: S/N of the Hi linewidth is low, but the 
shape of the line-profiles are very similar. 

PGC 16274: Nangay observations had a poor base-line, 
making a comparison between Ha and Hi line- 
widths difficult. 

PGC 19767: Ha data exhibits a strong asymmetry in 
the line profile that is not observed in the H 1 data. 

PGC 23598: Ha and Hi line- widths are equivalent and 
both profiles are asymmetric. The detailed shapes 
of the profiles, however, are markedly different. 
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Fig. 7. — Left Panel: Narrowband Ha image of PGC 38268 with overlay of the DensePak fiber locations in which Ha was detected. 
Aside from the Ha flux missed at large radii due to the limited field of view of Densepak, the IFU successfully detected Ha across the 
galaxy even at locations in which the imaging data does not yield discrete Ha detections. Center Panel: The DensePak Ha velocity field 
overlayed on the narrowband Ha image of PGC 38268. Right Panel: Normalized flux density of the PGC 38268 Ha line profile. Jagged 
dashed line is the raw line profile obtained by convolving the velocity field model and the Ha flux from the narrowband image. The heavy 
solid line is this same profile convolved with a Gaussian with 20 km s _1 width (which corresponds to the measured mean emission line 
width before instrumental correction was applied). The thin solid line is the DensePak line profile obtained by summing the individual 
Gaussian fits. The simulated profile has virtually the same width as the data and shows a qualitatively similar asymmetry profile. 



PGC 23913: Neither line profile shows evidence of sig- 
nificant rotation. 

PGC 24788: H I data has a slightly broader line width 
and is double-peaked while the Ha profile is not. 

PGC 26140: There is a strong mismatch between Ha 
and Hi line profiles, with little Hi signal at the 
source redshift, but a possible, narrow peak blue- 
shifted by 300-350 km s _1 . Within the Nancay 
beam lies PGC 26140 (2'A separation), an Ellip- 
tical with redshift (V=8570±33 km s" 1 ; LEDA), 
size, and magnitude similar to those the face-on 
spiral (V=8760±141 km s _1 ). Though our central 
Hi velocity (8575 km s _1 ) is closer to the opti- 
cal velocity of the Elliptical, we expected it to be 
gas-poor and our 21cm line detection should not 
be confused by it. PGC 26140 was not detected 
at Arccibo (Giovanardi & Salpeter 1985) with a 
quoted 3ct upper limit of 1.6 Jy km s _1 for a 300 
km s _1 wide flat-topped profile, considerable lower 
than our detected integrated line flux of 3.1 Jy km 
s _1 . We can conclude that the Hi mass for this 
galaxy lies far below the predicted mass, ~ 10.3 log 
Mhi, based on its luminosity and type (Roberts & 
Haynes 1994). 

PGC 28310: Ha and Hi line- widths are equivalent and 
both profiles are asymmetric. The detailed shapes 
of the profiles, however, are markedly different. 

PGC 31159: Hi and Ha line-profiles are equivalent, 
but the Hi profile is double-peaked while the Ha 
is not. 

PGC 32091: No Hi flux was detected. 



PGC 32638: Nangay observations have a poor base- 
line, making a comparison between Ha and H I line- 
widths difficult. 

PGC 33465: Nangay observations have a poor base- 
line, making a comparison between Ha and H I line- 
widths difficult. 

PGC 36925: No Hi flux from the source was detected. 

PGC 38908: Hi and Ha line-widths are very similar. 
The only notable difference is the slightly broader 
width of the Ha profile. 

PGC 39728: Hi and Ha line-widths are very similar. 
The H I profile displays a strong asymmetry that is 
not observed in the Ha profile. 

PGC 46767: Hi profile has a low S/N, but the Ha 
and H I profiles and profile-asymmetries are similar 
nonetheless. W 2u measured from Nangay profile is 
significantly larger than one would expect from the 
profile due to the baseline noise. 

PGC 55750: The spatial coverage of the DensePak 
pointings were insufficient to cover the disk of the 
galaxy (Observations were performed for only 2 of 
3 planned pointings). 

PGC 56010: Neither line-profiles show much evidence 
for rotation. The spatial coverage of the DensePak 
pointings were insufficient to cover the disk of the 
galaxy (Observations were performed for only 2 of 
3 planned pointings). 
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PGC 57931: Nangay observations have a poor base- 
line, making a comparison between Ha and H I line- 
widths difficult. 

PGC 58410: Hi and Ha line-profiles are very similar. 

PGC 72453: Ha and Hi line- widths are equivalent and 
both profiles arc asymmetric. The detailed shapes 
of the profiles, however, are markedly different. 

6. SUMMARY AND DISCUSSION 

We used the Principal Galaxy Catalog to choose a sam- 
ple of 39 disk galaxies which span a range in luminosity, 
surface brightness and type. For this sample, we gath- 
ered Ha emission-line data for 39 nearly face-on spiral 
galaxies with the DensePak IFU on WIYN. These obser- 
vations provide a fairly uniform spatial sampling of the 
ionized gas out to radii of approximately 3 scale lengths 
with a spectral resolution sufficient to yield emission-line 
velocities with a centroiding accuracy of ~ 2 km s _1 . 

For one of these galaxies, PGC 38268, we collected 
a deep narrow-band Ha image. By comparing the 
DensePak data to the narrow-band image, we found that 
the 2.8" DensePak fibers and 3x longer exposures en- 
abled us to detect discrete line-emission over a much 
larger filling factor than apparent in the narrow-band 
images. Despite the different detection thresholds, we 
used the PGC 38268 narrow-band image convolved with 
a model of the velocity field obtained from the DensePak 
data to construct a spatially-integrated, Ha line-profile 
which accounts for all Ha gas in the galaxy. A line pro- 
file constructed using only the DensePak data, despite 
not accounting for all the ionized gas, still was able to 
produce a comparable, integrated line-width. Qualita- 
tively, the two Ha emission-linc-profilcs exhibited the 
same width and asymmetry. Having demonstrated that 
DensePak observations can be used to produce a real- 
istic line profile for PGC 38268, we generated spatially 
integrated Ha line-widths for all galaxies in our sample. 

We have also measured the spatially-resolved velocity 
dispersion for ionized gas in our sample. For the 6000 
IFU spectra which were well-fit by a single Gaussian (and 
therefore less sensitive to beam smearing), we remark 
that under the assumption that the thermal line-width 
for ionized gas at 10 4 K is a = 9.1 km s _1 (Osterbrock 
1989), and the intrinsic line- width is a = 3 km s _1 , then 
following the quadrature formalism of Jimenez- Vicente 
et al. (1999), the turbulent motions for Hn in our sam- 
ple ranges from 7.7 to 24.2 km s _1 (80%). This range 
may seem large, but is in keeping with the amplitude 
and range seen in Fabry-Perot studies of much smaller 
samples of galaxies (cf. Jimenez- Vicente et al. 1999). 
Perhaps surprisingly, significant variations and similar 
amplitudes are seen in the molecular and neutral gas - 
albeit in small samples. For instance, Combes & Bec- 
quaert (1997) find 6 and 8.5 km s _1 dispersions for CO 
in NGC 628 and NGC 3938, respectively. While van der 
Kruit & Shostak (1982, 1984) have remarked on the con- 
stancy of H I velocity dispersions in disks ranging from 
7 to 10 km s _1 , this is based on radio-synthesis obser- 
vations of three galaxies. One of these galaxies, NGC 
1058, was re-observed by Dickey et al. (1990) at higher 
spectral resolution. They found a lower H I velocity dis- 
persion of 5.7 km s _1 . Both Shostak & van der Kruit 



(1984) and Kamphuis & Sancisi (1993) have also noted 
H I velocity dispersions in two late-type, nearly face-on 
disk decrease from 10-13 km/s in the inner disk (what 
Kamphuis & Sancisi call the "turbulent regions of the 
optical disks") to 6-8 km/s in the outer disk. While the 
physical mechanisms driving the velocity dispersions of 
the different ISM phases may indeed be different (cf. Mc- 
Kee & Ostriker 1977, and Jog & Ostriker 1988), our pri- 
mary points here is the variation of turbulent motions 
from galaxy to galaxy and within galaxies in all phases 
may be more significant than previously supposed. We 
will return to this issue in the context of the ionized gas 
in latter papers in this series. 

We gathered available, spatially-integrated Hi line- 
profiles from the literature and augmented them with 
new 21-cm observations from the Nangay radio telescope. 
We have presented new Hi line-profiles for 23 galaxies, 
along with measures of their line- widths and Hi fluxes. 
These data allow us to compare line-profiles of ionized 
and neutral hydrogen. For galaxies with sufficient S/N 
in their Hi profiles, we found excellent agreement be- 
tween the width of the profiles indicating that we are 
sampling the flat portion of the velocity fields in Ha for 
most of this sample. However, the sample of H I and Ha 
line-profiles showed significantly different profile shapes 
indicating that we are observing different distributions 
in the ionized versus neutral gas. 

We have shown that integral field spectroscopy does an 
admirable job of obtaining integrated line-widths. Line- 
profiles from Nangay have S/N ratios far less than the 
S/N ratios for line-profiles constructed from DensePak 
data. Yet the Nangay radio telescope has a collecting 
area more than 800 times greater than WIYN, and we 
used an average on-source integration two times longer 
per galaxy in Hi than for Ha. Even if we ignore the 
different signal-to-noise ratios of the lines, and account 
for the fact that radio telescopes are roughly 100 times 
less expensive per m 2 than optical telescopes (of compa- 
rable steerability; van Belle et al. 2004), and that Hi 
observations can be carried out during daylight hours, 
the enhanced cost-effectiveness of DensePak-to-Nangay 
observations is greater than a factor of 5. Since our data 
was taken, Nangay has undergone an extensive upgrade 
which has increased the sensitivity by a factor of ~ 5, 
which brings the cost effectiveness roughly to unity. High 
signal-to-noise ratio H I line profile observations using the 
Green Bank 43m and Arecibo 305m radio telescopes have 
comparable cost effectiveness to our DensePak observa- 
tions (e.g. Haynes et al. 1998; Haynes et al. 1999). 
Yet the true forte of Ha integral field spectroscopy is 
the high spatial resolution of resolved velocity fields and 
the simultaneous observations of multiple atomic species. 
Typical radio-synthesis H I maps require 10 to 20 hours 
of integration on Westerbork or the VLA, with most ra- 
dio observations using beam sizes of 15 arcsec or greater. 
This is characteristically five times longer integrations at 
5 times lower resolution than DensePak. H I and Ha ob- 
servations probe different ions and physical scales, and so 
remain complementary. However, with the relatively re- 
cent advent of bi-dimensional optical spectroscopy with 
DensePak and related IFUs, it appears that efficient kine- 
matic studies within the inner (optically bright) regions 
of galaxies are now squarely in the optical regime. 
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Fig. 8. — Leftmost Column: Images of the 39 sample galaxies with overlays of the DenscPak fibers in which Ho is detected. Central 
Column: Ha velocity profiles from the DenscPak observations. The Ha flux densities are plotted against heliocentric velocities. No units 
are given on the flux densities because nights were generally clear but not photometric, nor were flux-standards observed. The smoother 
profiles arc sums of Gaussian fits to individual fiber Ha profiles. These "Sum of Gaussian" profiles agree well with the raw data. Flux 
from fibers with centers that fall within a fiber radii are averaged. W20, the width of the Ha lines at 20% the peak flux is illustrated 
by a horizontal line in the plots. Rightmost Column: The Hi flux densities from Nangay 21cm observations (where available) versus 
heliocentric velocities. The smooth curve is again the Ha "Sum of Gaussian" profile scaled to the Hi peak flux. W20, the width of the 
Hi lines at 20% the peak flux is illustrated by a horizontal line in the plots. 



